Pseudogap in Doped Mott Insulators is the Near-neighbour Analogue of the Mott Gap 
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We show that the strong coupling physics inherent to the insulating Mott state in 2D leads to a 
jump in the chemical potential upon doping and the emergence of a pseudogap in the single particle 
spectrum below a characteristic temperature. The pseudogap arises because any singly-occupied 
site not immediately neighbouring a hole experiences a maximum energy barrier for transport equal 
to t 2 /U, where t is the nearest-neighbour hopping integral and U the on-site repulsion. The resul- 
tant pseudogap cannot vanish before each lattice site, on average, has at least one hole as a near 
neighbour. The ubiquituity of this effect in all doped Mott insulators suggests that the pseudogap 
in the cuprates has a simple origin. 



In a Mott insulator with one orbital per site, each unit 
cell contains an odd number of particles but the Fermi 
energy lies in the middle of a gap. In contrast, band 
insulators contain an even number of electrons per unit 
cell and the Fermi energy lies atop a full band. Con- 
sequently, Mott and band insulators are not adiabati- 
cally connected. Nonetheless, most Mott insulators tend 
to order antiferromagnetically below some temperature, 
Tjv. As a consequence, it is standard to view the Mott 
state simply as an antiferromagnet in which the unit cell 
has doubled. On this view, the insulating properties of a 
Mott insulator are equivalent to those of a band insulator. 
Band insulators possess rigid bands and hence doping 
only creates quasiparticles at the Fermi level. That this 
picture fails fundamentally in the parent cuprates, which 
are all antiferromagnetic Mott insulators, is immediately 
evident from optical conductivity experimentsQ,0 which 
reveal that even for T ^> Tjv, a gap of order 2eV ex- 
ists and doping leads to a massive reshuffling of spectral 
weight from 2eV to the Fermi energy. These experiments 
lay plain that what is missing in the antiferromagnetic 
reduction is Mottness itself: 1) in the absence of mag- 
netic ordering (T > Tjv), a charge gap exists in the sin- 
gle particle spectrum, 2) each electronic state in the first 
Brillouin zone has spectral weight both above and below 
the charge gap, and 3) the sum rule that each single- 
particle state carries unit weight is satisfied^ only when 
the spectral function is integrated across the charge gap 
not simply up to the chemical potential as in a band in- 
sulator. A consequence of (3) is that in the Mott state, 
the traditional notion that the chemical potential demar- 
cates the boundary between zero and unit occupancy fails 
fundamentally. This failure is central to Mottness. 

The breakdown of the traditional band insulator sum 
rule in the cuprates is well described[3j by the Hubbard 
model in which the on-site energy for double occupancy 
leads to a charge gap at half-filling. Such local on-site 



physics dominates the insulating behavour at half filling. 
In the lightly doped regime, S ss 0, effective interactions 
of longer range come into play as neighbouring sites now 
become correlated. If on-site correlation leads to a charge 
gap at half-filling, it is certainly a possibility that nearest- 
neighbour correlations for 8 ~ 0, for example, might lead 
to a suppression of the density of states at the chemical 
potential as well. In fact, it is well documented that 
all the underdoped cuprates possess a pseudog ap0 m 
the single-particle spectrum. However, the origin of this 
phenomenon is unknown. As the pseudogap does not 
appear to be a true T = phase and the pseudogap 
line joins continuously to the Mott insulator, proposals 
which require broken symmetry [fj 0.0 are difficult to 
reconcile with the Mott state. 

Without any assumption as to the nature of the ground 
state, we show that the the electron spectral function 
for the 2D Hubbard model contains a dip at the Fermi 
energy which results in a pseudogap in the single par- 
ticle density of states (DOS). The pseudogap remains 
pinned at the Fermi level in the underdoped regime but 
moves above it at an intermediate doping level, as is seen 
experimentally^. The pseudogap is fundamentally tied 
to local correlations on neighbouring sites much the way 
Mott gap arises from on-site physics. 

The starting point for our analysis is the Hubbard 
model with nearest-neighbour hopping matrix element 
t and on-site Coulomb repulsion U. We base our strong 
coupling analysis on a two-component composite basis 
ip with Vio-(i) = £i<7 and ip2a-(i) = Via- an d its as- 
sociated Green function S(i,j,t,t') — {{^ia]^] a )) = 

9(t - H)({$ia{t),il>) a (f)}), where = c* CT (l - ?w) 
and rjia- = Ci a ni^ a . Here, Ci a = r\i a + annihilates 
an electron on site i and rii is the number operator for 
site i. The basis ip exactly diagonalizes the on-site in- 
teraction and hence serves as a natural starting point 
for a strong-coupling analysis. To overcome the stan- 
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dard truncation problems inherent in the use of Hub- 
bard operators, we adopt the following procedure. First, 
project all new operators that arise from the Heisen- 
berg equations of motion of the Hubbard operators onto 
the Hubbard basis. Second, write the self energy exactly 
in terms of the remaining operators which are now or- 
thogonal to the Hubbard basis. Third, use local methods 
in the spirit of dynamical mean- field theory (DMFT)0 
to calculate the resultant electron self energy. To go be- 
yond the single-site treatment indicative of dynamical 
mean-field theories^fji we adopt the two-site expansion 
proposed by Mancini and Matsumoto^J m which the 
self energy is determined self-consistently from a two-site 
Hubbard cluster embedded self-consistently in an inter- 
acting bath. As all orientations of the two-sites are con- 
sidered, the electron spectral function will be momentum 
dependent. Self-consistent cluster methods which are ex- 
act as the limit of an infinite cluster appear to be rapidly 
convergent, providing accurate results for the thermo- 
dynamics of the ID and 2D half-filled bands[9j and in 
fact constitute the accepted methodology^] for treat- 
ing strongly correlated systems. Hence, an implementa- 
tion of the Hubbard operators coupled with DMFT-type 
technology places the limitations not on truncation in 
the equations of motion but on the accuracy of the im- 
purity solver and the size of the finite cluster. As the 
complete procedure is detailed elsewhere |9Hl2|. we men- 
tion only that in contrast to the work of Matsumoto and 
Mancini [TH. we required that for a fixed filling in the lat- 
tice, the chemical potential of the cluster equal that of 
the lattice. 

Using this procedureplH^. we report first the doping 
dependence of the chemical potential. Two distinct pos- 
sibilities arise: 1) the chemical potential remains pinned 
and mid-gap states are generated or 2) the chemical po- 
tential jumps across the Mott gap. Our results shown in 
Fig. (JJJ demonstrate that the chemical potential jumps 
upon hole or electron doping, indicating an absence of 
mid-gap states. The Magnitude of the jump is set by 
the Mott gap which is fully developed at T = 0. Even 
for U — 4t, the inset on the right shows that the chem- 
ical potential resides in the LHB for n = 0.95. While 
at some finite temperature, the chemical potential may 
appear to evolve smoothly, A/i ^ as the doping in- 
creases and hence no mid-gap states are present. Exact 
results in the ID Hubbard model as well as Quantum 
Monte-Carlo simulations 14] in 2D also reveal a chemi- 
cal potential jump upon doping and hence no mid-gap 
states. However, unlike ID and 2D, in d — oo, the chem- 
ical potential remains pinned^(| upon doping as mid- 
gap states emerge. A chemical potential jump requires 
a large imaginary part of the self-energy at the chemical 
potential. From the inset in Fig. (JJJ, we find that 
is initially large in the underdoped regime and acquires 
the characteristic u> 2 dependence in the overdoped regime 
indicative of a Fermi liquid. Consequently, the method 
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FIG. 1: Doping dependence of the chemical potential in the 
2D Hubbard model computed using the local cluster approach 
for T = 0.15t (dashed line) and T = 0.07i (solid line). The 
inset on the left shows the imaginary part of the self energy 
evaluated at a Fermi momentum (0.3,2.10) for n = 0.97, 
(0.3, 1.84) for n = 0.8 and (0.3, 1.06) for n = 0.3, whereas the 
inset on the right contains the density of states for n = 0.95 
for U = At, U = 8t, and U = 12t. 



we use here is capable of recovering Fermi liquid the- 
ory in the overdoped regime. Experimentally, whether 
Afi vanishes or not appears to be cuprate dependent. 
For example, in La2- :E Sr ;E Cu04 (LSCO), the che mi- 
cal potential remains pinned roughly at OAeV above the 
top of the LHB, while for Nd 2 _ x Ce ;E CuO ± i (NDCO)[lf| 
and Bi 2 Sr 2 Ca 1 _ ;E R x Cu208+j / (BSCO) [TIE Hill, the 
chemical potential jumps upon doping and scales roughly 
as S 2 as obtained here. Because stripes require A^i = 0, 
they have been invoked [2 lj to explain the origin of the 
mid-gap states in LSCO. The pseudogap in the under- 
doped cuprates has also been attributed^J to stripes. 
However, because A/i ^ for both NDCO and BSCO, a 
requirement for stripe formation, if the pseudogap has a 
universal origin in the cuprates, stripes are not its cause. 

To address the origin of the pseudogap, we focus on 
the doping dependence of the electron spectral function, 
-S(Sn + 2Si2 + S 2 2)/k), shown in Fig. ©. Six features 
are evident: 1) the chemical potential moves further into 
the LHB as the filling decreases, 2) no coherent peaks 
exist near the chemical potential for n — .97, 3) each 
state in the FBZ has spectral weight both above and be- 
low the chemical potential as dictated by Mottness, 4) 
the Mott gap remains intact but moves to higher energy 
as the doping increases, 5) at (n, 7r), the UHB carries 
most of the spectral weight regardless of the filling, and 
6) a dip exists in the spectral function at the chemical 
potential for n — 0.97 but is absent for n — 0.60. In 
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the underdoped regime, the characteristic width of each 
k state is of order t and even much larger near (71", 0). 
Such broad spectral features in the underdoped regime 
are seen experimentally^^ and arise in this context be- 
cause ^ as shown in Fig. JIJ. As a consequence, 
there is no sharp criterion for unit occupancy of each 
state in the FBZ. Because the spectral weight at each 
momentum is unity, however, and each state lives both 
below and above the chemical potential, the charge car- 
ried by the piece of the state lying below the chemical 
potential is less than unity. 

Is the dip in the spectral function shown in Fig. 
for n — 0.97 real? Displayed in Fig. J3J) is the DOS for 
T = 0.25t and T = O.Cm for several fillings. As is ev- 
ident, no local minimum of DOS exists at the chemical 
potential at high temperature, T = 0.25t. Features which 
emerge even at high temperature are the reshuffling of 
spectral weight from above the charge gap to below as 
the filling is changed and also a movement of the Mott 
gap to higher energies. Note that even at n — 0.30 the 
Mott gap is still present, though almost all of the spectral 
weight now resides in the LHB which closely resembles 
the non-interacting density of states. This is further ev- 
idence that we correctly recover Fermi liquid theory as 
71 — > 0. At low temperature, the lower panel of Fig. 
© demonstrates that a pseudogap forms in the DOS for 
S « 0. The vertical line at indicates that the pseudo- 
gap occurs precisely at the chemical potential. Similar 
qualitative results based on a cluster method have been 
obtained by Maier, et. al.[22^. except their pseudogap is 
slightly displaced above Ep. In contrast, in the analy- 



FIG. 3: Density of single particle states for T = 0.25t and 
T = 0.07t, U — 8t for the fillings shown. No pseudogap 
exists at high temperature. At low T, a pseudogap emerges 
and remains pinned at the Fermi level but moves above at 
an intermediate doping level. In the overdoped regime, the 
pseudogap vanishes entirely and a non-interacting system is 
recovered. 



sis of Haule, et. al. 23], the DOS has a negative slope 
through Ep (as dictated by the proximity to the Mott 
gap) but never acquires a local minimum at Ep indica- 
tive of a true pseudogap. Because the pseudogap exists 
below some characteristic temperature and vanishes at 
higher doping, the result obtained here is non-trivial and 
highly reminiscent of experimental trends Q. What is its 
origin? The inset in Fig. indicates that for a fixed 
filling, the pseudogap vanishes as U increases and scales 
as t 2 /U. This suggests that the pseudogap is tied to 
short-range correlations and hence explains why it is ab- 
sent in d = 00 01- While J e g oc 0{t 2 /U) corresponds to 
the energy scale for antiferromagnetic spin fluctuations, 
such fluctuations cannot inhibit hole transport. In fact, 
Maier, et.al.j^ have shown that even if antiferromag- 
netism is killed, the pseudogap still persists. Further, we 
have found that J e g is only weakly doping dependent for 
< x < .25 and in fact vanishes at x ~ .8. Hence, the 
resolution of the pseudogap problem lies elsewhere. The 
energy scale, t 2 /U also describes any transport process 
in which the intermediate state is doubly occupied. Such 
processes are captured by our approach as a result of the 



4 



coupling of the two-site cluster to the interacting bath. 
Consider placing a single hole in a Mott insulator. Un- 
like a site neighbouring the hole, a singly-occupied site 
two lattice sites away must temporarily doubly occupy 
one of its neighbours if it is to move to the hole. The 
energy for this two-step process is t 2 /U. Sites further 
away experience an energy barrier with a higher power 
of t 2 /U. Hence, t 2 /U is the largest energy barrier for 
hole transport once a Mott insulator is doped. Because 
some sites experience no energy barrier, the single par- 
ticle density of states exhibits only a suppression, a sig- 
nal that hole transport involves virtual excitations to the 
UHB. This pseudogap cannot vanish before each site has 
on average one hole as its immediate neighbour, roughly 
x = .25 for a square lattice. Hence, the pseudogap is of 
the form t 2 /UP(x), where P{x) determines the proba- 
bility that hole transport involves double occupancy and 
consequently, is a steadily decreasing function of x. 

Additionally, it is precisely two-step (or three-site) 
hopping that makes the single-particle low-energy spec- 
tral weight increase faster [5^] than 2x. To showthat we 
recover this result, we compute the high and low spectral 
weight by integrating the DOS from the energy which 
minimizes the DOS to oo (— oo for electron doping) and 
from p, to that fixed energy, respectively. The results 
shown in Fig. Q (which have been normalized per spin) 
demonstrate that the initial spectral weight in the UHB 
which is 1/2 at n = 1 all moves to low energies as the 
filling decreases as is observed experimentally[l|,|3|- The 
same is true for electron doping (n > 1). Further, the cur- 
vature of the low energy spectral weight (LESW) is pos- 
itive indicating that the LESW grows faster[3|, |25j than 
2x. The growth in excess of 2x arises entirely from vir- 
tual excitations between the LHB and UHB and points to 
an iseparability of the low and high energy scales. Such 
behavior is absent from a band insulator (see Wni in 
Fig. That Mottness leads to such a drastic devi- 

ation from the non-interacting result is a direct conse- 
quence of each state having spectral weight both above 
and below the chemical potential (see Fig. J2J). Our 
finding that the three-site terms lead to an inseparability 
of low and high-energy scales resonates with the recent 
work of Kirkpatrick and Belitz|26| who have shown that 
three-body terms are ubiqutious in strongly correlated 
electron systems and lead to breakdown of a true low- 
energy description. 

Without global symmetry breaking, spin-charge sep- 
aration, or pairing, we have shown that lightly doped 
Mott insulators posess a pseudogap which arises entirely 
from nearest-neighbour correlations. The pseudogap is 
ubituitous because any singly-occupied site not immedi- 
ately neighbouring a hole experiences an energy gap for 
transport equal to t 2 /U. This generic phenomenon of- 
fers a simple resolution of the pseudogap problem in the 
cuprates. 

This work was funded by the Petroleum Research Fund 
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FIG. 4: High (Wh) and low (Wl) spectral weight as a 
function of filling. Wni is the spectral weight in the non- 
interacting system. 
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